Coherent Time Evolution and Boundary Conditions of Two-Photon Quantum Walks 
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Multi-photon quantum walks in integrated optics are an attractive controlled quantum system, 
that can mimic less readily accessible quantum systems and exhibit behaviour that cannot in general 
be accurately replicated by classical light without an exponential overhead in resources. The ability 
to observe time evolution of such systems is important for characterising multi-particle quantum 
dynamics — notably this includes the effects of boundary conditions for walks in spaces of finite size. 
Here we demonstrate the coherent evolution of quantum walks of two indistinguishable photons using 
planar arrays of 21 evanescently coupled waveguides fabricated in silicon oxynitride technology. We 
compare three time evolutions, that follow closely a model assuming unitary evolution, corresponding 
to three different lengths of the array — in each case we observe quantum interference features that 
violate classical predictions. The longest array includes reflecting boundary conditions. 



Random walks describe stochastic motion of a particle 
around a discrete space and are widely applied as a st a _ 
tistical tool in areas ranging from genetics to economics. 
The quantum mechanical version — quantum walks — is 
an interference phenomenon of quantum particles that 
exhibits distinctly different dynamics. Quantum walks on 
many graphs exhibit ballistic propagation of the walker's 
probability distribution, while many trials of a classical 
random walk will build a distribution localised around 
the initial position pQ — a 1-dimensional quantum walk 
is a typical example. For certain tasks, quantum walk 
dynamics are beneficial, and have inspired new quan- 
tum algorithms [2] and approaches to universal quan- 
tum computing by realising quantum walks on a finite 
set of graphs [3] ; in this setting, quantum walks are tra- 
ditionally thought of as being simulated on a quantum 
computer. Quantum walks are also a physical phenom- 
ena in their own right, describing for example evolution 
in solid state systems [4j [5] and transport phenomena in 
biomolecules [6 — all such systems are finite and there- 
fore exhibit either periodic, absorbing, reflecting or cyclic 
boundary conditions which have been treated theoreti- 
cally, for example [7J [8] . 

There are two main models of quantum walks [1] — 
continuous time and discrete time, of which single walker 
examples have been demonstrated experimentally using 
trapped ions [9HTT]: optical resonators [12]; NMR [HlfHj : 
single photons in bulk p~5] [16] and fibre optics [T7J [18] ; 
and laser light in coupled waveguide arrays [19]. The 
observable dynamics of single walkers can be described 
by either single particle quantum mechanics or classical 
wave mechanics [20], as demonstrated by optical quan- 
tum walks using laser light. When multiple indistin- 
guishable particles undergo a quantum walk, correlations 
occur that cannot in general be explained or mimicked 
with classical physics, without either sacrificing visibility 
in observable features [21 — quantified by an inequality 



relating correlated outcomes — or introducing a factori- 
ally increasing number of experiments [22 or resources. 
The latter approach has been recently adopted with op- 
tical fibres and laser light [23]. Two particle quantum 
walks [21] have been demonstrated with photons guided 
in integrated optics: Indistinguishable photons have lead 
to quantum correlations in a planar [24] and a circular 
[25] array of evanescently coupled waveguides; polarisa- 
tion entangled quantum walks have recently been demon- 
strated and shown to emulate the quantum interference 
statistics of fermions, bosons and the intermediate regime 
[H[27]. 

Here we measure three distinct time steps of a two- 
photon continuous-time quantum walk in three arrays 
of 21 waveguides, with identically designed coupling and 
propagation parameters, integrated on a single chip. The 
length of each array z is varied, to provide the mechanism 
for observing multi-photon dynamics in time [z = c-t for 
speed of light c in the structure). The longest time step 
exhibits boundary effects from photons being reflected at 
the outer waveguides. In all three cases, we observe quan- 
tum interference that violates an inequality that relating 
correlations between different waveguides and assumes 
classical light as the input [21]. 

Photons propagating in a planar array of N evanes- 
cently coupled single mode waveguides realise a quan- 
tum walk on a one dimensional lattice (see Fig. [I] (a) for 
N=21). Assuming nearest-neighbour coupling from the 
evanescent fields between waveguides j and j db 1 (am- 
plitude ftj,j±i = ttj±i,j) an d a waveguide propagation 
constant f3j, the Hamiltonian is modeled with 

N 

H = c^2(3 j a]a j + K di i _iat_ 1 a j + Kjj+ia] +1 aj, (1) 

where a\ and a 7 are the bosonic creation and annihilation 
operators for a photon in waveguide j and for uniform 
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arrays Kij and j3j are constant. 

Time evolution of a given state \ijj(t)) is determined 
by the unitary operator U(t) = exp(-iHt) according to 
\ip(t)) = U(t) \i/)(to = 0)). The transition amplitude from 
the initial state of a particle in waveg uide j (a] |0) = |1) ■) 
at time to = to the state of a particle in waveg- 
uides j' (at, |0) = at time t is then given by 
Ujfj(t) = (l| •, exp(— iHt) yielding the single par- 
ticle transition probability 



(2) 



which also describes the photon density distribution one 
expects in experiments with coherent light. 

Quantum correlations occur in the output state of the 
quantum walk of multiple indistinguishable particles [21] . 
From injecting a photon pair into waveguides k and j, the 
probability to detect one photon in waveguide j' coinci- 
dent with detecting one photon in waveguide k' is given 
by the two-photon correlation function [28] 



rj$(*) 



i 



1 + Sj^k' 



\Uj>At)Uk>,k(t) + u k ,j(t)Uk,j>(t)\ 2 . 

(3) 

Non-classical features emerge from interference of the 
complex probability amplitudes inside the | • | 2 . Ex- 
pression (|3| with unitary operator U is the theoretical 
description to which we compare our experimental re- 
sults. The correlation function of distinguishable par- 
ticles is treated using classical probability theory, multi- 
plying and summing the possible single photon transition 
probabilities: V'^) {t) = Pj' ,j(t)p k > , k (t) +Pj' 1 k(i)pk'j(t). 
To distinguish between two-photon quantum interference 
experiments and any classical treatment of light, refer- 
ence [21 defined the inequality 



V Jt k{t) = 



rk-(*)n, fc (*)-r;, fc (t)<o 



(4) 



which is only violated by the interference of indistinguish- 
able photons [24] . 

The quantum walk Hamiltonian ([I]) is implemented 
with waveguide arrays fabricated in silicon oxynitride 
(SiO^N^) with a designed refractive index contrast of 
4.4%. This high index contrast enables fabrication of mi- 
cron sized single mode waveguides which are pitched at 
2.8/im within the coupling region in order to achieve suf- 
ficient mode overlap for nearest neighbour coupling. The 
waveguides are designed with constant height of 0.6/mi 
and width of 1.8/xm. Three waveguide arrays with cou- 
pling region lengths z — 350/im, 700/im and 1050 jam 
were fabricated on a compact 5 mm long chip (Fig. [TJ 
and measured to simulate three values of time evolution 
parameter t = z/c. 

Identical, horizontally polarised 808nm photons gener- 
ated from type-I downconversion source (see Appendix) 
are launched into each of the three arrays using a butt- 
coupled array of polarisation maintaining, single mode 
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FIG. 1: Waveguide structure for measuring time evolution. 

(a) Graph representing the planar array of 21 waveguides. 

(b) Chip with the three measured waveguide arrays. One can 
see the three input waveguides for each array, followed by the 
narrow coupling region and the spreading to the output at the 
top of the chip, (b) Schematic of the waveguide structure; the 
waveguides run parallel within the coupling region and bend 
to the input and output separation uniformly. 



optical fibre with a standard pitch of 250 /im; this ad- 
dresses the input of the array by bending the input 
waveguides to a separation of 250 jam (with a maximum 
bend radius of 600/im). The output of the array is 
similarly addressed by bending the output waveguides 
to a separation of 125 jim and butt-coupling the out- 
put chip to a second 250 jam pitch array allowing corre- 
lated two-photon detection with commercially available 
fibre-coupled avalanche photo-diode single photon count- 
ing modules, across even and odd labeled waveguides. 
Two-photon correlations were measured by recording 
pairs of detection within a 5ns window as a coincidence 
event. Coincidence logic was performed with field pro- 
grammable gate arrays. Photon-number resolving de- 
tection was performed non-deterministically using fibre 
beam splitters. 

The single photon distribution Pj>j(t) was charac- 
terised with 810nm laser light injected into the central 
waveguide (input waveguide 11). Fig. [2] shows the sin- 
gle photon distribution measured with a power intensity 
meter, for the three different array lengths in compari- 
son with the theoretical prediction calculated according 
to equation (|2|. 

Two photon correlations Tjr^{t) were measured from 
injecting indistinguishable photon pairs into the two pairs 
of input waveguides ({11, 12} and {10, 12}) that pro- 
vide the most overlap throughout the device, for each of 
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FIG. 2: Single photon evolution, (a) Simulated probability 
distribution for a single photon walk starting in the middle 
of the array at waveguide number 11. (b-d) Three measured 
intensity distributions / normalized to the total intensity for 
the arrays of (b) 350 /im (c) 700/xm and (d) 1050 /im. 



the three simulated time evolutions. Indistinguishability 
of the photons was verified by performing a Hong-Ou- 
Mandel type experiment [29 with an average visibility 
of V = 98.4 ± 0.5%. The correlation matrices for inputs 
{10, 12} for the three propagation lengths are shown in 
Fig. [3] (a)-(c) (correlations for inputs {11, 12} given in 
the Appendix). 

Each of the 121 measured two-photon coincidence pos- 
sibilities (from all 231 possible combinations) were mea- 
sured over 15 minute integration periods and are cor- 
rected for relative detector efficiency. We quantify the 
agreement between the theoretical calculations r*v k ,(t) 

and the experimental data T^f p k ,(t) with the similarity 

2 

\^3 f ,k f Y L j',k' ' L j',k') I l^j' ,k' 1 3' ,k' 2^j> ,k' 1 3' ,k" 

In order to derive the values of the entries of the Hamil- 
tonian ([!]) of our waveguide structure we use standard 
coupled mode theory [30 according to which the prop- 
agation constant of waveguide j is given by /3j = no- 
where n is the effective refractive index and A the vacu- 
um wavelength. The coupling constant n between two 
identical waveguides can be deduced from the propaga- 
tion constants of the two normal modes f3 e and f3 Q of the 
coupled system as n = (/3 Q — /3 e )/2. We use a waveguide 
mode solver (FIMMWAVE) to calculate the correspond- 
ing propagation constants. The uniform array with the 
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Ztot 


S, input 11,12 


S, input 10,12 


350/xm 


425/xm 


94.0 ± 0.2% 


90.0±0.2% 


700/im 


775/mi 


90.3±0.1% 


90.1±0.1% 


1050/mi 


1125/xm 


89.6±0.1% 


90.7=b0.1% 



TABLE I: Similarities for different time steps and 2-photon 
input combinations. 



designed waveguide dimensions can be modeled with cou- 
pling and propagation constants of k = 0.00489/im _1 and 
f3 = 11.40462/im -1 . The spreading regions at the input 
and output add an effective coupling length of 75/im to 
each array length, so that the total coupling lengths are 
given by z to t = Uot * c = 425/im, 775/im and 1125/im. 

Maximising the similarity of the measured two-photon 
distribution with theoretical predictions we obtain a 
Hamiltonian (Eq. ([T])) with 



Pi 



11.385 for i = 1,21 
11.403 for i = 2,20 
11.397 else, 



0.0049 for i = 1,20 
0.0058 for i = 2, 19 
0.00513 else. 



where m = = 

This suggests fabrication deviation since boundary ef- 
fects like an additional self coupling term added to the 
propagation constant of outermost waveguides are a mag- 
nitude smaller than the derived perturbations [31 . The 
resulting similarities between theory and experiment are 
summarised in Table |U 

We verify the quantum behaviour of the photons at all 
three time steps through the violation of the classical in- 
equality Q as depicted in Fig. |3jg)-(i). The array length 
of 700 jam shows the strongest violation with a maximal 
value of 128 standard deviations, while the longest ar- 
ray (1050 /im) shows comparatively small violations of 
a maximum of 45 standard deviations. This behaviour 
is in agreement with a theoretical study of the time de- 
pendence of the violation (Fig. [4|. We calculate the 
violations for the theoretical correlation matrices (for all 
121 measurable two-photon coincidences) and plot the 
maximum value. 

Within the time interval relevant for our experiment 
the maximum violation decreases with time as the pho- 
tons distribute over the waveguides. For longer times, 
as the photons reached the boundaries of the array and 
are back-reflected, the maximum violation shows an os- 
cillatory pattern. The local minima and maxima for the 
maximum violation reflect the changing phase relation 
between the matrix elements of the unitary time evolu- 
tion operator Uj t k. 

In summary, we have observed the time evolution of a 
two-photon continuous time quantum walk on an array 
of 21 evanescent ly coupled waveguides. The three mea- 
sured time-steps agree with theoretical prediction made 
assuming coherent unitary evolution; this is evidence of 



4 



(a) | (b) | (c) | 




FIG. 3: Two-photon data. Experimental [(a)-(c)] and theoretical correlations [(d)-(f)] for detecting simultaneously one photon 
in waveguide j f and one in waveguide k' when injecting photon pairs in waveguide 10 and 12. (g)-(i) shows the violations of 
equation Q in units of standard deviations. 



V 

0.04 
0.03 
0.02 
0.01 




425 775 1125 



2500 



tc 



FIG. 4: Time evolution of maximum violation V m ax for in- 
put 10 and 12 in comparison with measured maximum vio- 
lations (points). Errorbars on experimental data are smaller 
than data points. Theoretical and experimental correlation 
matrices used to calculate the violations are normalized as 
probability distributions. 



coherent evolution of multi-photon states within inte- 
grated optical circuits. However, full characterisation 
of time evolution can be achieved via state tomogra- 



phy and as the size of quantum walks increase, schemes 
that approximately determine the quantum process us- 
ing resources efficiently could be implemented [32]. We 
have observed that the magnitude with which two-photon 
quantum interference violates classical prediction varies 
over the three measured time steps. Our numerical sim- 
ulations support these fluctuations, and specifically pre- 
dict that it is possible to have zero violations, despite 
genuine two-photon interference. In further array de- 
signs the pitch and width of the outer waveguides can be 
adapted to model various boundary gradients. 
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Appendix 

Experimental results for input waveguides 11 

and 12. Here we present the experimental results for a 
second input combination when injecting photon pairs in 
neighbouring waveguides numbered 11 and 12 (see Fig. 1 
(c) in main text). Fig. [5] shows the correlation matrices 
and violations measured. Quantum interference results 
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(a) | (b) | (c) | 




FIG. 5: Experimental [(a)-(c)] and theoretical correlations [(d)-(f)] for detecting simultaneously one photon in waveguide f 
and one in waveguide k' when injecting photon pairs in waveguide 11 and 12. (g)-(i) shows the violations of equation (4) in 
the main text in units of standard deviations. 



in distinctly different correlations compared to the input 
combination of waveguides 10 and 12. In all three time 
evolutions we observe a generalized bunching behaviour, 
the two particles are likely to go to the same side of the 
array. Fig. [6] compares theoretically predicted time evo- 
lution of the maximum violation with our measurements. 
Compared to the input {10, 12}, the photons wavefunc- 
tions start overlapping at an earlier time, so violations 
can be observed at an earlier time (around t-c = 130/im). 
Despite genuine two-photon interference, the maximum 
violation goes down to zero at t • c = 230/im. 

Two-photon down conversion source. The pho- 
ton pairs were produced via type I spontaneous paramet- 
ric down conversion in a 2-mm-thick y 2 nonlinear bis- 
muth borate BiE>306 crystal. A continuous- wave diode 
laser (404nm) pumped the crystal with 60 mW to create 
degenerate photons at a wavelength of A = 808nm. In- 
terference filters (3.1nm) ensured spectral indistinguisha- 
bility before we launched the light into two polarization- 




FIG. 6: Time evolution of maximum violation Vmax for in- 
put 11 and 12 in comparison with measured maximum vio- 
lations (points). Errorbars on experimental data are smaller 
than data points. Theoretical and experimental correlation 
matrices used to calculate the violations are normalized as 
probability distributions. 



maintaining fibres (PMF), which were butt-coupled to 




FIG. 7: Schematic of the type-I down conversion source used 
to generate (unent angled) photon pairs. 



